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Abstract 

On the basis of a moment method, general solutions of a linearized Boltzmann equation for a 
normal Fermi system are investigated. In particular, we study the sound velocities and damping 
rates as functions of the temperature and the coupling constant. In the extreme limits of col- 
lisionless and hydrodynamic regimes, eigenfrequency of sound mode obtained from the moment 
equations reproduces the well-known results of zero sound and first sound. In addition, the mo- 
ment method can describe crossover between those extreme limits at finite temperatures. Solutions 
of the moment equations also involve a thermal diffusion mode. From solutions of these equations, 
we discuss excitation spectra corresponding to the particle-hole continuum as well as collective 
excitations. We also discuss a collective mode in a weak coupling case. 



PACS numbers: 52.35.Dm, 67.25.dt, 67.85.Lm 



I. INTRODUCTION 



In discussing collective excitations of quantum degenerate gases, there arises a distinction 
between hydrodynamic and collisionless regimes. Let lo be a frequency of a collective mode 
and t be a mean-collision time. The hydrodynamic regime is characterized by lot <C 1, 
while the collisionless regime is characterized by lot ^> 1. The present paper is devoted to 
detail analyses of collective excitations in both limiting regimes and in the crossover regime. 
We deal with a sound propagation in a population balanced normal Fermi gas. 

In general, the viscosity calculated in the hydrodynamic regime diverges at the absolute 
zero temperature (T = 0), and hence it was considered that sound could not propagate 
at T = 0. In 1957, using the Fermi liquid theory Landau predicted that a sound 
propagation could occur in liquid 3 He even at very low temperatures owing to the mean- 
field interaction [2|. This new type of sound was called zero sound, which differs from first 
sound propagating because of a small dissipation achieved by local equilibrium. 

Khalatnikov and Abrikosov conducted detail calculations of the dispersion relation and 
derived the sound attenuation coefficient of the zero and first sound modes [3J]. Abel et. 

n 

al. confirmed the existence of zero sound in the liquid 3 He [4|]. They also observed a 
crossover between the zero and first sound modes, and measured temperature dependence 
of the sound velocity and of the sound attenuation coefficient. The theoretical prediction 
based on the Landau's Fermi liquid theory [3j agreed with the experimental data. Until now, 
investigations of liquid 3 He have been conducted in detail, and the results are summarized 
in many text books [5]. The Landau's Fermi liquid theory is also summarized in standard 
textbooks [6|. 

With the realization of the Bose-Einstein condensate as a turning point [7], vigorous 
studies of ultracold atomic gases have been conducted. Ultracold atomic gases have flexibil- 
ities such as controllability of an interaction parameter using the Feshbach resonance. These 
systems open new windows to investigate phenomena that were difficult and impossible to 
study in the liquid helium and superconductors. 

As mentioned below, a number of studies on the collective modes have been also reported 
in ultracold Fermi gases. Dipole oscillations were studied in collisionless and collisional 
regimes in two component 40 K gases s), 0]. Experiments of the collective excitation in 
the BCS-BEC crossover regime have been performed using 6 Li gases [lO, H, 12, 13]. In 
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particular, the reference [13j reported the sound velocity in the BCS-BEC crossover regime. 

Collective excitations in Fermi gases have been also investigated theoretically. The dipole 
mode Jjj] and the quadrupole mode 15] were analyzed making use of the moment method. 



Zero sound with arbitrary spin was investigated 



in trapped gases extensively and intensively 
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6]. Bruun et. al. studied collective modes 



18l . Il9|. Tosi 's g roup studied collective 



modes by solving the Boltzmann equation numerically 2CJ, |2l|, |22j, 123J. They studied the 



dipole mode o 
2(1 



uations 



trap [22 



two component trapped g clSGS clS db function of a collision rate in some sit- 



2l| . and the crossover between zero and first sound modes in the cigar-shaped 



23| . Recently, collective excitations in the unitarity limit were studied 



As noted earlier, after a publication of the path-breaking work by Landau 
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26]. 



2JL a classic 



paper 
sound 



3V Khalatnikov and Abrikosov studied the crossover between zero sound and first 
3j. They treated a mean-collision time as temperature-dependent, but approximated 
other quantities by those at T = 0. When the system is in the collisional regime at finite 
temperatures, a sound velocity within this analysis reaches that of first sound evaluated at 
T = 0. This simple approximation is appropriate as long as we restrict discussion to the 
Landau's Fermi liquid theory, since this theory focuses on quasiparticles at sufficiently low 
temperatures. In the experiments in liquid 3 He [4], however, the temperature dependence 
of first sound velocity has been observed, although it was small. In usual atomic gases, 
furthermore, first sound has a significant temperature dependence, and hence the simple 
treatment mentioned above is invalid at high temperatures. It is therefore necessary to 
discuss the crossover between the zero and first sound modes with a more efficient treatment 
valid at wide range of temperatures. 

Brooker and Sykes attempted to analyze the general solution of the linearized Boltzmann 



equation to investigate the crossover of sound propagation [27]] . They expanded a deviation 
from local equilibrium in terms of spherical harmonic functions, and introduced different 
relaxation times for different spherical harmonics. They, however, used an approximation 
only appropriate in the low temperature regime, and hence the analysis was not appropriate 
at high temperatures. They, furthermore, could not obtain the explicit solution because of 
the computational difficulty at that time, although they gave an equation to be solved. 

Although the physics on the crossover between the zero and first sound modes has been 
understood to some extent, some issues shown above still remain. With the developments 
of recent experimental techniques in ultracold atomic gases and of theoretical methods, it is 



3 



meaningful to revisit the study of sound mode in Fermi systems with a modern approach. 
In the present paper, we investigate the crossover between zero sound and first sound over 
wide parameter ranges with a single theoretical framework. For this purpose, we analyze 
the general solution of the linearized Boltzmann equation using the moment method. 

The moment method is suitable for describing the crossover of collective excitations be- 
tween collisionless and hydrodynamic regimes. So far in cold atomic gases, the moment 
method has been used to study characteristic collective excitations in trapped systems, such 



as monopole, dipole, quadrupole and scissors modes [28|, |29|, |3Cj, |3l|. With a use of same 



tec hnique, collective modes in atomic gases with internal degrees of freedom were also stud- 



ied 



32 



33( | . The moment equation for the uniform system, however, has not been solved to 



study sound modes in quantum many-body systems. This is one of the new points in the 
present paper. 

In the Landau's Fermi liquid theory, central players are quasiparticle, and hence the- 
oretical studies have been done in the very low temperature regime. Those studies are 
based on the Landau Boltzmann equation for quasiparticles. In dilute quantum gases, on 
the other hand, real atoms are central players. These systems are described by the Boltz- 
mann equation for real particles. Although the Boltzmann equation has the same form as 
the Landau Boltzmann equation, this equation is applicable up to the high-temperature 
Maxwell-Boltzmann gas regime. Studies of sound mode in ultracold atomic gases thus need 
a single theory which can be applied up to high temperatures. The present study can solve 
those issues. 

The contents of the present paper are summarized as follows: (a) The spectrum of the 
sound mode obtained from the general solution shows the crossover between zero sound and 
first sound. In collisionless and collisional regimes, the sound velocities reproduce the results 
calculated in each limiting regime. This method also offers the frequency of first sound with 
the temperature dependence. This result cannot be obtained by a standard approach such 
as given in Ref. 3j. (b) The results of the moment method include a thermal diffusion mode, 
(c) The moment method reproduces not only a collective mode of the sound propagation, 
but also the particle-hole continuum, (d) In a weak coupling case, the crossover between 
zero sound and first sound has a different character from the crossover in a strong coupling 
case. 

The present paper is organized as follows. Section [III deals with one of the main topics 
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of the present paper. Making use of the moment method, we will derive moment equations. 
In Sec. IHIt we calculate the sound velocity and damping rate of first sound in the hydro- 
dynamic regime. Section IIVI gives detailed analysis of zero sound in the low temperature 
regime. In Sec. |V] we will show results of moment equations, and analyze the crossover 
between zero sound and first sound. Section [VI] is devoted to discussion. Section IVIII gives 
summary and conclusion. We devote Appendices [A] and [B] to derive relaxation times. In 
Appendix [A] we calculate transport coefficients in the hydrodynamic regime based on the 
Chapman- Enskog method. In Appendix [B], we evaluate relaxation times associated with the 
transport coefficients. We compare the mean collision rate with these relaxation rates: the 
viscous relaxation rate and the thermal conductivity relaxation rate. Appendix IC1 describes 
a standard analysis of the random-phase approximation. 



II. LINEARIZED BOLTZMANN EQUATION AND MOMENT EQUATION 

In this paper, we consider two component atomic Fermi gas interacting with s-wave 
scattering. We assume a population balanced gas of two spin components with the same mass 
m. The equation of motion for distribution functions within a semiclassical approximation 
is described by the following Boltzmann equation: 

dUV d ^ t] + £ ■ V r / CT (p,r,t) - V r U a (r,t) ■ V p / CT (p,r,i) = Z con [f a ], (1) 

where an index o = {t? 1} represents spin component. U a (r,t) is the contribution of a 
mean-field interaction given by Ucr(r,t) = gn^^r, t). n a (r,t) represents the local density. 
The interaction strength g is given by g = 4-n-h 2 a/m, and a is the s-wave scattering length. 
The collision integral T co \\[fo\ on the right hand side of Eq. ([I]) is given by 

x [1 - [1 - f_ a (2)} /- CT (3)/ CT (4) - / CT (1)/_ CT (2) [1 - /_(3)] [1 - / CT (4)] , 

(2) 

where f a (i) is f a (i) = f a {pi,r,t). 

We shall linearize the distribution function around static equilibrium (denoted by 
/°(p,r)), using f a {p,r,t) = /°(p, r) + 5/ CT (p, r, t) = /°(p, r) + 6f a (p, r, t) + 8f {j>, r, t). 
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Here, 5/ CT (p, r, £) is the local equilibrium distribution function linearized around static equi- 
librium /cr(p, r, t) = /°(p, r) + Sf a (p, r, t), and Sf'(p,r,t) denotes departure from local 
equilibrium. 

The local equilibrium distribution is determined by the condition 1 co \\ = 0, and is given 

by 

r, t) = p — - — , , (3) 

exp < (3(r, t) [p — mv(r, t)] /2m >z~ x {y, t) + 1 

where the local fugacity z a (r, t) is 

z a (r,t) = exp{/3(r,£) \fi a {r,t) - gn„ a (r, *)]}. (4) 

/9(r, t) is a local temperature 9(r, t) = fc B T(r, £) = l//3(r, t), ?v(r, £) is a local density, /x CT (r, £) 
is a local chemical potential, and v(r, £) is local velocity. These local variables depend on 
position and time. 

It is convenient to write fluctuations of the distribution function around static equilibrium 

as 

df° 

*/«r(P,M) = -Q^V„{p,T,t). (5) 



The factor f CT (p, r, t) is the averaged extra energy of particles around equilibrium 34|. Writ- 
ing 5f a (p,r,t) and 5/£(p,r,t) as 

df° df° 
5f a (p,r,t) = ^P a (p,r,t), 6f (p,T,t) = ^5u a (p,r,t), (6) 

one also has ^(p, r, £) = £ CT (p, r, £) + Su a (p, r, £). 

Using Eq. (jSJ) in the Boltzmann equation ([Tj), one obtains the following equation: 



di/ ff (p,r,t) p p 
H Vz/ ff (p, r, t) - Vgdn^fr, t) 



Xcoll[/ CT ]. (7) 



at m m 

We shall apply a relaxation time approximation to the collision integral. The relaxation time 
t is a characteristic time with which a system reach local equilibrium. In this approximation, 
the collision integral can be reduced to 

ZconL/y = = (8) 

We now linearize the local equilibrium quantities as 9(r, t) = 9° + 59(r, £), fi a (r, t) = n° a + 
6ji a (r, t), and v(r, t) — v° + 5v(r, £) = Sv(r, t), where 9° = k B T°, chemical potential /jP, and 



velocity v° = represent static equilibrium. The linearized local equilibrium distribution 
function is then given by 



where 



P a (p,r,t) = a a (r,t) + b(r,i) • p + c(r,t)p 2 



a a (r,t) = -[p gn _ a -p n a ]59(r,t)-5fi a (r,t) + g5n^ a (r,t), 

b(r,t) = Sv(r,t), 

c(r,t) = -^6e(r,t). 
2m 



(9) 

(10) 

(11) 
(12) 



We now look for the plane wave solution of the linearized Boltzmann equation, represent- 
ing as 



z/ CT (p, r, t) = u a (p, q, u;)e* (q ' r wt \ 5n a (r, t) = 5n a (q, to)e 



where resultant functions v a and 8v a are also written as 



i(q-r— u}t) 



u a (p, q, to) = U a (p, q, to) + 6v a (p, q, u), 

z> CT (p, q, lo) = a ff (q, lo) + b(q, lo) ■ p + c(q, uj)p 2 . 



(13) 

(14) 
(15) 



The linearized Boltzmann equation with the relaxation time approximation is thus given by 
df° a 



del 



p-q\ , x . pq x ( s 

lo ) v a {p, q, lo) + —— gdn^ a {q, to) 



m 



1 df° 



Here, the density fluctuation 5n a (q, to) is 



f dp df° 

5n a (q,u) = a a {q,uj)W a>0 + c(q,uj)W (T! 2 + J ^^ ^^'^ ^ 



where we define 



(27rft) 3 de° a ' 

From Eqs. f|T4l) - (|T71) . the linearized Boltzmann equation is reduced to 



dfl 
del 



lo 



p ■ q 



m 



[a a (q, lo) + b(q, lo) ■ p + c(q, ^)p 2 + Su a (q, 



LO 



, pq 

+ 9 

m 



dp df° 

a^ a (q,Lo)W^ }0 + c(q,u)W^ 2 + / J^Kf de° " 6u - a ^ ^ 



(16) 



(18) 



.1C 

t oe" 

(19) 



We now discuss the general solution of the linearized Boltzmann equation, a main topic 
of the present paper. As derived in the above, the linearized Boltzmann equation with the 
averaged extra energy around the static equilibrium v a {p) is reduced to 

[(" - — ) ^(P) + —9Sn J = \v a - {a. + b ■ p + cp*)] . (20) 

We do not write explicitly q, and u in v a and 5n CT , since these dependences are not important 
for further calculation. We shall use the viscous relaxation time given in Eq. (1B9I) (derived 
in Appendix [A] and [Bj as the relaxation time r, because the density fluctuation is the most 
strongly coupled to the viscous relaxation. This approximation to the collision integral is 
a good one in the vicinity of local equilibrium. In the collisionless regime lot ~^s> 1, this 
approximation can describe zero sound, because the collision term can be neglected owing 
to the large value of r. We remark that the relaxation time evaluated by a small correction 
from static and local equilibrium could be quantitatively different from the actual relaxation 
time in the collisionless regime lot ^> 1 and also in the intermediate regime lot 1. 
We expand the fluctuation in terms of the spherical harmonics as 

oo I 

^(P) = EE K n A P )Pr(cos9)e m t (21) 

1=0 m=-l 

Multiplying Eq. (!20l by e~ tm '^ and integrating it over 0, we have the following linearized 
Boltzmann equation: 

df® \( pq 



1=0 



(22) 



C^T (cos 9) - (a a + b • p + cp 2 ) 5 mfi 

.1=0 

One finds that only the mode with m = is coupled to the mean-field potential. This is 
due to the isotropic interaction. This mode (m = 0) corresponds to the longitudinal wave. 
In anisotropic interactions, there also exists the mode with m/fl, such as transverse zero 
sound with m — 1. Since we consider the crossover from the longitudinal zero sound to 
the longitudinal first sound, we only take the mode with m = 0. Let us use the notations 
^ =0 (p) = v a,n{p), and P z m=0 (cos6 l ) = Pz(cos#), for simplicity. It is also useful to define the 
following moment: 

<"""*.'> s / T0W% pn "° M - (23) 
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The density fluctuation is expressed as Sn a = (v a ,o)- 

Multiplying Eq. (1221) by p n Py (cos 9) and integrating over 9 and p, we obtain the moment 
equation given by 

21 — lm 21 + 6m m 

= ~ -(P n v«,i) + - (civW^n + cW^n+z) 5;,o + -bW a , n+1 5i A . (24) 
r r r 

We have made use of an orthogonality relation 

/ d9sm9P l {cos9)P v {cos9) = 7^T7 S W> ( 25 ) 

and a recurrence formula for the Legendre polynomials 

cos0P*(cos0) = i±l^ +1 (cos0) + ^pY^_i(cos0). (26) 

The moments associated with p°Po(cos9), pPi(cos9) and p 2 Po(cos9) correspond to num- 
ber of particles, momentum, and the energy, respectively. The collision integral vanishes 
when we take these moments, because of the conservation law. Equations determining co- 
efficients a a , b, and c are then given by 

Ko) - a a W^ - cW a , 2 = 0, (27) 
E(K,i>-^,2) =0, (28) 

^((p\o>-8A, 2 -*,4) =0. (29) 

We used an assumption that directions of the velocity and of the sound propagation are 
parallel b || k, where b is related to the velocity through Eq. ffTTj) . As a result, one obtains 
coefficients given by 

a ° = W^ o) " e >> + e ^ w^ M > (30) 
6= <p y t,i) + (pn,i) ni N 
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where = ^2 (W^ — W^ 2 /W a fi). Finally, we obtain the following moment equation: 

I q 



21 -lm 



T 

i 

T 



:<P%i> 



w a>n 1 fw a , n w* 2 w a>n+2 w a , 2 



e 



(T.O 



fw a , n w a , 2 w^, 2 W a , n+2 W^, 2 



+ 
+ 



r6 
i 1 

i 



(v<7,o)8i,o 

(V-a,o)8l,0 



W a)n+ 2 ~ 

w a>n+1 



w a , 



W- a , 

((p 2 ^T,o) + (p 2 n,o)) 8i,o 



(33) 



T W^ 2 + W ij2 

One can obtain the eigenmode by solving this eigenvalue problem; however, equations 
are not closed even if higher moments are taken into account. We shall truncate an equation 
at sufficiently high moment, which does not affect the spectrum of the collective mode of 



interest. Note that this equation is not the same one derived in Ref. 
is much simpler than that in Ref. 



271 ] . The equation 



(T33D 



271 ] . We do not use 



many relaxation times as in Ref. 



23, 



271 ] added an extra equation to make 



but a single relaxation time is introduced. Reference 
a closed set of equations, but we do not need an extra equation. In next two sections, we 
grasp sound velocities and damping rates in the two limiting regimes: hydrodynamic and 
collisionless regimes. 



III. FIRST SOUND 



We solve the linearized Boltzmann equation in Eq. f)19p in the hydrodynamic regime in 
the present section. When we take the zeroth, first, and second moments of the Boltzmann 
equation, the collision integral vanishes because of conservation laws. As shown in Eqs. 
( 1A46|) and (1A47I) in Appendix [Aj one can obtain a closed set of hydrodynamic equations 
including dissipative terms. The hydrodynamic equations in terms of the moments are 
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written as 

uj [a CT (q, uj) W a>0 + c(q, W„, 2 ] = b(q ^ ' V ^, (34) 
> <^ wb q,w ^ 2 - a a (q,uj) q-c q,w q 

a v 

[a_ CT (q, c^WU.o + c(q, u;)W_ CT , 2 ] q 1 - *4 W 2 b(q, w) = 0, (35) 
w [a T (q, w)W^ Ti2 + a l (q,uj)W l<2 + c(q,a;)(W Ti 4 + W^)] 

- b ^ q :^' q (ty T ,4 + ^4,4) - 4mm 2 c(q, cu)Tg 2 = 0. (36) 

Relations 59(q,uj) = fce5T(q, uj) and b(q, a;) = — v(q, uj) are used, and the velocity v(r, t) 
is assumed to be parallel to a vector q. 

The above equations can be written in terms of physical quantities: the density Sn a (q, uj), 
the velocity v(q, uj) and the energy 5E(q_,uj), whose quantities are given by 

5n a (q, uj) =a CT (q, uj)W afl + c(q, uj)W^ 2 , (37) 

v(q,w) = - b(q,w), (38) 

<5£(q, =o T (q, w)W T , 2 + fli (q, w)W i)2 + c(q, w)(W T)2 + W i>2 ). (39) 

Density fluctuations of an in-phase mode <5n to t(q, w) = 5nj (q, a;) + <foj(q, a;) and of an out- 
of-phase mode <m_(q, a;) = Sn^(q, uj) — 5nj_(q, a;) exist because of the two component system. 
The hydrodynamic equations in terms of these quantities are written as 

= uj5n_{q,uj), (40) 
2W 2 

= uj5n tot (q,uj) + - — q-v(q,w), (41) 
6m 

uj-i—) q . v(q, W ) + 2-^^(q, ") " ^^tot(q, a;), (42) 

0= (cj-4mm 2 Tg 2 i?)5E(q,w) + ^q-v(q,cu)-4mm 2 Tg 2 J R— <5n tot (q,cj), (43) 

3m lVo 

where i? is defined as 1/R = 2(W4 — W^/Wo), an d the assumption of the population balanced 
gas W n = = W 7 ^ is used. We note that the out-of-phase mode is decoupled from the 
hydrodynamic mode composed of the total density, the velocity, and the energy, in the 
population balanced gas. 

Solving the secular equation for fluctuations (5n tot (q, a;), qv(q, a;), SE(q, uj)), one obtains 
an equation 

F 1 (uj) + F 2 (uj)=0, (44) 
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where 



F 1 (u) ^ 3 - u^y-m - gW*), (45) 
F 2 (u) = ( ^! + AiKm 2 Tq 2 R) + -mTq^R - gW 2 ) . (46) 



v w 2 ";3" \w 2 / 

We omit terms of second and higher order in transport coefficients k and rj. These coefficients 
are assumed to be small in the hydrodynamic regime. 

A frequency u of a collective excitation can be separated into a real part Q and an 
imaginary part T with being a damping rate: uj = Q — iT. Undamped solutions satisfy 
Fi(jY) — 0- Frequencies Q ^ obtained from Fi(Q) = are 



n^±n = J w >-J w ^. (47) 

1 3W 2 m K J 

A mode Q = which is a thermal diffusion mode also exists, and it will be discussed in Sec. 
VIII. 

In the weak coupling limit at T — > 0, the frequency is given by ~ qvp/\/3, where vp is 
the Fermi velocity given by vp = [h/ m)^^ 2 N tot /V) 1 ^ . N tot is the total number of particle, 
and V is a volume. In the strong coupling limit at T — ► 0, on the contrary, the frequency is 
given by ~ q^J g N to t/ (2mV) = qv-p a/<7Pf/3, where pp is the density of state at the Fermi 
energy given by pf = 2m(37r 2 N to t/V) 1 ^ 3 /(2nh) 2 . 

Damping rates and transport coefficients are assumed to be small in the hydrodynamic 
regime. The term F 2 (Q — iT) including transport coefficients can be then approximated by 
F 2 (Q), so that we reduce Eq. (jHJ) as Fi(Q — iT) + F 2 (Q) =0. As a result, one obtains the 
damping rate r± as 

. F 2 (n±) r/q 2 KTq 2 m 2 



n±(n±-n T ) w 2 w 4 - g w 2 

where we use Fi(u) = uj{uj — Q + )(u — fi_). 



(48) 
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IV. ZERO SOUND 



This section discusses a sound mode in the collisionless regime. We start with the lin- 
earized Boltzmann equation based on Eq. 



[("-^ cos 0)1/^ (cos 



1=0 



/I COS I 



— (-OOSflj^o) 

K,o> <^T,l>+<pn,l> pcos 



«=0 



(49) 



where the simplified notations ^ = °(p) = v a,n{p) an d -P™ =o (cos0) = P/(cos6>) are used. We 
assume conservation laws only for number of particles and for momentum as in Ref. 

Oct = (f<T,o)/W<T,o, & = Z)(P I/ ff,i)/X)^ ) 2, and c = 0. 
Equation (1491) can be reduced to 



1 0/° 1 



cos 6* 



i I COS I 



1=0 



cos 6* 



1 



ir del pq a 

a COS V 

m 



ITU) 



%Tpq I m 



x(^-<x,o> 

p cos 



(50) 



Multiplying this equation by 1 and p cos 6, and integrating over the momentum p, we obtain 
the following two equations: 



1 



iTq/m Wo 
1 A h2 



A . 



iTq/m W 2 



(pvi) = gAi, 2 + 



Ao: 



(51) 
(52) 



3 iTq/m W 2 J \ ' iTq/m W 

where W n = W^ n = W^ n . We consider the in-phase mode (jp n vi) = (p"z/y) + (p n ^i,i)- The 
coefficient A n j is defined as 



dp df° 

(2vr/i) 3 de° 



p" 



cos 1 6 



cos 9 



iTUJ — 1 



(53) 



iTpq/m 

where the spin index in df°/de° is omitted. 

Let us consider zero sound in the low temperature regime. We impose the temperature 
dependence only to the relaxation time, and simply evaluate the coefficient A n j as the value 
at T = 0. From Eqs. ( JBTl) and ( J52l) . the following dispersion relation can be obtained: 



1 -i— + W"(s) 
s 



gpF + i h i3s" (s — is") 



0. 



(54) 
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The function W°(s) is the Lindhard function given by 

8 + 1 



1 



2 



s-1 



(55) 



s and s" are defined as s = {ituj — l)/(irvpq), and s" = l/(rvpq), respectively. The 
dispersion relation (154"]) is first derived by Khalatnikov and Abrikosov [3]. The notations 
used here follows in Ref. [35 ]. 

The frequency of the collective mode in collisionless limit tuj ^> 1 is given by 



i + 9Pf i - 



2vpq 



In 



f2 — fp? 



0. 



(56) 



where we omit s" in the dispersion relation ( 1541) because it is small. This result can be also 
obtained by the random phase approximation discussed in Appendix O The frequency is 
given by uj = qv F {l + 2exp [— l/(gp F )]} in the weak coupling limit gpp <C 1; the frequency 



is given by uj = qvp^gpp/3 in the strong coupling limit gpp 1. 

The damping rate T = —lm{u) in the collisionless limit tuj ^> 1 is obtained by the 
following way. We expand the dispersion relation (|54|) to first order in s" and V since these 
are small. As a result, we obtain the damping rate Y as 



1 

T 



(gp F + l + 3sg)(sg-l) 
9Pf(9Pf + 1 - si) 



(57) 



where so = VL/v^q. This is seen in Ref. 35] . 



The frequency and the damping rate in the low frequency regime lot 1 can be also 
evaluated based on the dispersion relation (|54l) . The Lindhard function W°(s) is approxi- 
mated as W°(s) ~ — l/3s 2 — l/5s 4 in this regime. When we consider the dispersion relation 
( |54j) with the first order of the relaxation time r, we obtain 

2 1 A 

(58) 



— J = \i.9PF 
v F qJ 3 



1) iuor. 

' 15 



As a result, the frequency of the collective mode is given by f2 = qv-p^J{l + gpp)/3. One 
obtains f2 ~ qv-p a/s'Pf/S in the strong coupling case gp F ^> 1. This corresponds to the 
frequency of the first sound at T = in Eq. (|4"T1) in the strong coupling limit. The damping 
rate can be approximately evaluated as T ~ 2r(vpq) 2 /15 from Eq. ( |58l) . This damping rate is 
consistent with the damping rate of the first sound in Eq. ( HHl) . if we impose the temperature 
dependence in Eq. (HHj) only to the relaxation time, and evaluate other quantities in Eq. 
PB]) as the value at T = 0. The viscous term — rjq 2 /W 2 alone contributes this damping rate. 
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V. RESULTS : SOUND MODE FROM cjt«1TOwt»1 



This section presents the results obtained by solving the moment equation (|33j) . We focus 
on the crossover from zero sound to first sound. 

The collisionless regime lot ^> 1 and the collisional regime lot <C 1 can be realized by 
controlling the temperature T. In the high temperature regime, atoms are colliding with 
each other frequently, so that the hydrodynamic regime lot <C 1 is achieved. In the low 
temperature regime, on the contrary, the Pauli blocking makes the phase volume where the 
atoms are scattered restricted, and hence the collisionless regime lot ^> 1 is realized. 

The coupling constant a = gN tot / iVe-p) also plays a role determining the collisionless and 
collisional regimes, where e-p is the Fermi energy given by e-p = [^ 2 /(2m)](37r 2 A^ tot /l / ) 2 / 3 . 
The mean-field potential is proportional to the coupling constant a, while the collision 
integral (or the relaxation rate) is proportional to a 2 . The collisionless regime lot ^> 1 could 
be realized in the weak coupling regime. In the strong coupling regime, on the contrary, the 
collision term (or the relaxation rate) is dominant compared with the mean-field potential, 
so that one would be in the collisional regime. From these points of view, we study the 
sound mode from lot <C 1 to lot ^> 1 as a function of the temperature T and of the coupling 
constant a. 

In Fig. [Tj, eigenvalue lo of the collective excitation is plotted as a function of T. We show 
a strong coupling case a = 15. We chose the wavenumber q = 0.05k-p. We take moments up 
to Z = 30 and n = 30 in this calculation, although less moments, for example up to I = 10 
and n = 10, reproduces the same result. 

Figured] (a) shows the phase velocity c defined by c = Q/q where fl = Re(u;). Figured] 
(b) shows damping rate T given by T = — Im(u;). Solid lines in Figs. [1] (a) and (b) represent 
the phase velocity and the damping rate obtained from the moment equation (1331) . Dashed 
lines in Figs. [1] (a) and (b) represent those of zero sound obtained from Eq. fl56l) and given 
in Eq. (1571) . respectively. Dotted lines in Figs. [1] (a) and (b) represent those of first sound 
given in Eq. ( 1471) and Eq. ( ]48|) . Solutions of the moment method coincide with asymptotic 
solutions in two limiting regimes: collisionless and hydrodynamic regimes. Note that the 
moment equations show the crossover between the zero and first sound modes as well as the 
temperature dependence of first sound. Corresponding behavior of our result is also seen 
in the experimental result in Ref. which reported temperature dependence of the sound 
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k B T/e F k B T/s F 



FIG. 1: Frequency and damping rate of collective excitation as a function of temperature. Figure 
(a) shows the phase velocity. Figure (b) shows the damping rate. Solid lines correspond to an 
eigenmode obtained from the moment equation. Dashed lines in (a) and (b) represent the phase 
velocity and the damping rate of zero sound, respectively. Dotted lines in (a) and (b) represent 
those of first sound. The coupling constant a = 15 is used. 

velocity and the amplitude attenuation coefficient of liquid 3 He. 

In the collisional hydrodynamic regime ujt <C 1, the dispersion relation in Eq. (15^|) first 
derived by Khalatnikov and Abrikosov cannot reproduce our results correctly, because 
the temperature dependence is imposed only to the relaxation rate (see also Eq. ( J58l) ). 
In Eq. ( H71) . W&j which is proportional to the pressure, strongly depends on temperature, 
and this brings temperature dependence of the sound velocity of first sound. (W 2 , which 
is proportional to the density, does not have the temperature dependence under a fixed 
volume.) 

As for the damping rate, the dispersion relation in Eq. (|58|) does not involve the con- 
tribution of the thermal conductivity. Even if we neglect the second term in Eq. ( 1481) . the 
dispersion relation in Eq. (1581) still does not reproduce our result, although the difference 
is quite small. The difference also comes from the temperature dependence of the pressure 
in the term W4, which is not involved in Eq. fl5^|) . It is, nevertheless, remarkable that the 
dispersion relation first derived by Khalatnikov and Abrikosov 3j excellently grasps the 
sound velocities and damping rates in both collisionless and hydrodynamic regimes. 

In turn, we plot the phase velocity c and the damping rate T of the collective mode as a 
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FIG. 2: Phase velocity in (a) and damping rate in (b) are shown as a function of the coupling 
constant a = gN tot /V£-p, fixing the temperature at T = 0.025£f- Solid lines in (a) and (b) show 
the results obtained by the moment method. The dashed lines in (a) and (b) are the phase velocity 
and the damping rate of zero sound, respectively. The dotted lines in (a) and (b) represent those 
of first sound, respectively. 

function of the coupling constant a in Figs. [2] (a) and (b). We show the low temperature 
case k^T = 0.025£f- Again, we choose the wavenumber q = 0.05/cf, and take moments up 
to / = 30, and n = 30. Solid lines in Figs. [2] (a) and (b) show the phase velocity and the 
damping rate obtained from the moment equation fl33l) . The dashed lines and the dotted 
lines in Figs. [2] (a) and (b) show the velocity and damping of zero sound (given in Eqs. (|56|) 
and (Iq"T|) ) and first sound (given in Eqs. flTTl) and fHS]) ). respectively. The crossover from 
zero sound to first sound can be clearly seen in this figure. 

From Fig. [2] (a), one can see that the phase velocity of zero sound is close to that of first 
sound in the strong coupling regime. This is due to the fact that the phase velocity c of 
zero sound is given by the same formula as that of first sound c ~ vp a/^Pf/S in the strong 
coupling limit at T = 0. Note that the mechanisms of sound propagation are quite different 
in two regimes. 

One could change the coupling constant a = gN tot / (Vep) by controlling a density N tot /V, 
or an interaction strength g through the Feshbach resonance. The Fermi energy is also a 
function of the density, i.e., ep oc (N tot /V) 2 ^ 3 , and hence the coupling constant a has a 
density dependence: a oc (N tot /VY^ 3 . 
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VI. DISCUSSION 



In this section, we discuss physical implication of results obtained from the moment 
equation. First, we discuss a hydrodynamic mode other than the sound mode. Second, 
excitation spectrum of the particle-hole continuum obtained from the moment equation is 
discussed. Third, we discuss the sound mode in a weak coupling case is made. Forth, other 
issues and future problems are discussed. 

A. Thermal Diffusion Mode 

In discussing the collective mode in the hydrodynamic regime, there usually exist five 
modes, corresponding to the particle number, the velocity and the energy. Two modes are 
the first sound modes ±f2 — iT, related to the particle number and the velocity of a certain 
direction, discussed in Sec. IIHI Other two modes are shear modes related to the velocity 
of remaining two directions. The other is the thermal diffusion mode T K related to the 
energy. Note that the shear modes and the thermal diffusion mode are purely damping 
modes. 

In the present paper, we assume that vectors b and q are parallel each other as treated 
in Sec. [Til arid Sec. IIHI This means that velocity of the fluid v is assumed to be parallel to 
the wavenumber vector of the collective mode q, and hence two shear modes are neglected. 
In this subsection, we discuss the thermal diffusion mode. 

In Sec. IHfl we noted that a collective mode Q = exists. Assume that the result is written 
as uj = —iT, and consider the term up to the first order in damping rate and transport 
coefficients in Eqs. f|4"5l) and (pEBj) . Setting Fi(-iT) + F 2 (0) = 0, we obtain damping rate of 
the thermal diffusion mode 

2 K Tm*q 2 W*(l-gWo) 

(w 4 - gw$)(w 4 w - w%y { ' 

Since the moment method provides the general solution of the linearized Boltzmann 
equation, the thermal diffusion mode should be also included. In Fig. [31 the damping 
rates corresponding to the thermal diffusion mode are plotted. In this calculation, we take 
moments up to I — 30 and n = 30, and chose the wavenumber q = 0.05fcp. The solid lines 
are the damping rates obtained from the moment equation (1331) . The dotted lines are the 
damping rates given in Eq. f[59l) . In Fig. 4 (a), the damping rate of the thermal diffusion 
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mode versus temperature is shown for the coupling constant a = 15. In Fig. 4 (b), the 
damping rate versus the coupling constant a is shown for the temperature k&T = 0.025£f- 
Parameters in Figs, (a) and (b) are the same as in Fig. [T] and Fig. [2J respectively. We 
confirm that the present moment method provides the thermal diffusion mode. 
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FIG. 3: Thermal diffusion modes are plotted as a function of temperature in (a), and of the coupling 
constant in (b). Solid lines in (a) and (b) are the results obtained by the moment method. The 
dotted lines represent thermal diffusion modes given in Eq. (|59p . 



B. Particle-Hole Continuum 

In discussing zero sound, one often uses the random phase approximation. The usual 
random phase approximation (see Appendix [C]) gives excitation spectra in the particle-hole 
continuum as well as a collective mode. We discuss excitation spectra at T = obtained by 
solving the moment equation from this point of view. 

In Fig. HI we plot real part of eigenvalues obtained from the moment equation fl33l) as a 
function of the wavelength q at T = 0. The real part of these frequencies is symmetric with 
respect to the g-axis, so that we show only the region Q > in FiglU We take moments up to 
I = 30 and n = 30 in the numerical calculation. The coupling constant a = gN tot / '(Ve-p) = 5 
is used. From Fig. HJ one finds that real parts of eigenvalues in the moment equation 
also include excitation spectra corresponding to the particle-hole continuum as well as the 
collective excitation. The gradient of the edge of the particle-hole continuum excitation in 
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FIG. 4: Eigenvalues Q = Re(w) as a function of the wavelength q at T = 0, obtained from the 
moment equation. 



this figure is seen to be 2 in our dimensionless units, which corresponds to ft — vpq in the 
real physical units. 

In the usual random phase approximation, the denominator of the response function is 
given by uj + e p — £ p+q , and hence spectrum includes u = p ■ q/m + q 2 /{2m) |36|. This 
feature brings the phonon excitation uj oc q at the long- wavelength regime q/kp <C 1, and 
the parabolic excitation u oc q 2 at q/k-p ^> 1, where kp is the Fermi wavenumber. Solution 
of the semiclassical Boltzmann equation only involves the denominator uj — p ■ q/m, as seen 
in Eq (1C4jl . and hence our calculation can reproduce only the phonon regime: Q oc q. 

We presented discussion of the particle-hole continuum, but we remark some issues shown 
in EFD] 



C. On the Weak Coupling Case 

The phase velocity of zero sound is always larger than the Fermi velocity when g > 0. 
The phase velocity of first sound, however, could be less than the Fermi velocity in the 
weak coupling case and at low temperatures. In such cases, the spectrum of the collective 
excitation is not necessarily pushed up above the particle-hole continuum. We discuss the 
results in such a weak coupling case. 
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We calculate phase velocities c = Ke(u)/q as a function of the temperature in a = 1. We 
chose q = O.Olkp. In the calculation, we take moments up to I = 11 and n — 11. A reason 
of truncation at the moderate moments is that it allows us to clearly see transitions of the 
each eigenvalue. At T ~ 0, zero sound is seen as a separate eigenvalue where c ~ Vp. At 
finite temperatures, we confirmed that the eigenvalue of the collective excitation is buried 
in the particle-hole continuum. In this calculation, we also found that the spectra of zero 
sound and of first sound are not continuous in contrast with the strong coupling case. This 
result suggests that a collective mode in a weakly coupling system has a different feature 
from that in a strongly coupling system. It is unclear that how the collective mode behaves 
in the weakly coupling system in the crossover regime, even if more moments are taken. In 
a separate paper, we will calculate the dynamic structure factor of a normal Fermi system 
at finite temperatures, and discuss this problem \^\ . 

D. Remarks and Future Problems 

Before closing this section, we make some remarks on results of the present method and 
propose future problems. 

We first note some issues on the present method. Even if we set the relaxation rate 1/r 
to be zero, the coefficient matrix of the moment equation ( 12^1) is not symmetric, although 
the matrix elements are real. The eigenvalues are thus complex, in general. Those damping 
rates, namely, the imaginary parts of the resulting eigenvalues, range from —Vpq to Vpq 
at T = 0, and hence the present moment method could not reproduce the results of the 
random phase approximation perfectly. We note that those damping rates increase as the 
temperature T or the coupling constant a increases. In addition, there exists an additional 
purely damped mode, which is not included in the sound mode or the thermal diffusion 
mode. This mode does not belong to the complex eigenvalues discussed above either. The 
damping rate of this mode could be negative at certain temperatures and certain coupling 
constant. Albeit the present moment method involves the issues mentioned above, we insist 
that the present method offers very intriguing studies on the collective mode over a wide 
range of parameters. 

In turn, we shall discuss the future problem from the physical point of view. 

The excitation spectra in the weakly interacting system is complicated as discussed in 
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Sec. |VT] C. The spectrum of the collective excitation buried in the particle-hole continuum 
in the crossover regime. The collective mode and the single particle excitations are strongly 
related in this crossover regime, and hence the effect of Landau damping could be important. 
One issue is how the feature of the collective mode remains or disappears in this regime. 
In a separate article, we will study this problem by the dynamic structure factor 37J. The 
effect of the Landau damping would be seen in the peak width of the dynamic structure 
factor. 

It is also interesting to solve the equation derived by Brooker and Sykes 27]]. The 
equations in Ref. 27j involves additional equation in order to close the moment equation. 
Ref. J27I, in addition, introduces different relaxation times for different moments. Such a 
treatment is complicated compared with our formulation, and it is not obvious how the 
additional equation affects our result. 

The Landau's Fermi liquid theory focuses on the low temperature property, since this 
theory is based on an idea that a lifetime of quasi-particles are sufficiently long at very low 
temperatures. For this circumstance, the crossover between the zero and first sound modes 
has been studied theoretically only within the low temperature approximation. In ultracold 
Fermi gases, the real-particle picture is also important in both a classical gas regime and 
a weakly interacting Fermi system. Our formulation allows one to describe such a system. 
The temperature T, the density N tot /V and also the interaction strength g are controllable 
with recent techniques in ultracold atomic gases. We expect that behaviors of the collective 
mode shown in the present paper would be observed in the experiments of ultracold Fermi 
gases. 

We now comment on the application of the present work to a strongly interacting Fermi 
gas near the unitarity limit. At sufficiently high temperatures of the Maxwell-Boltzmann 
regime, real particles are important, and thus a gas is described by the Boltzmann equation 
with an energy dependent cross section 25j. In contrast, in the low temperature regime 
above the superfluid transition temperature, the system may be described by the Landau's 
Fermi liquid theory for quasiparticles. One expects a crossover from quasiparticle picture 
to real particle picture with increasing temperature, which may be studied by the moment 
method developed in the present paper. Explicit determination of the range of s-wave 
scattering length as well as of the temperature, where the kinetic equation analysis based on 
the long-living quasiparticle picture is valid, will require many-body calculation for a Fermi 
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gas near the unitarity limit 



38] 



VII. SUMMARY AND CONCLUSION 



The moment method is suitable for describing the collective mode from collisionless to 
collisional regimes with only a relaxation time approximation. We solved the linearized 
Boltzmann equation for a normal Fermi system using this method, and obtained the general 
solution. We discussed the crossover between the zero and first sound modes as a function of 
the temperature and the coupling constant. We found that an eigenfrequency of a collective 
mode obtained from the moment equations reproduces the sound velocity and the damping 
rate in the crossover regime as well as both collisionless and collisional limiting regimes. 

Through the analysis of the moment equation, we found that the moment method provides 
the thermal diffusion mode. We also discussed the excitation spectra of the particle-hole 
continuum, and the sound mode in a weak coupling case. We finally made remarks on the 
present method and proposed future problems. 
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APPENDIX A: CHAPMAN-ENSKOG METHOD AND TRANSPORT COEFFI- 
CIENTS 



In this section, we give a derivation of transport coefficients in a degenerate Fermi gas 
based on the Chapman-Enskog method. The result in this section will be used to evaluate 
the relaxation time in the next section. 



The transport coefficient in the Landau's Fermi liquid was first calculatec 
and Khalatnikov 



39l | . Afterwards, it was analyzed in several papers 
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jy Abrikosov 
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43). The 
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Chapman- Enskog method was first generalized to quantum gas es by Uehling and Uhlen- 



beck 
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451 ] . In this Appendix, the analysis is based on Ref. 461 ] . 



Following the standard procedure, we define the following hydrodynamic physical quan- 
tities: 



density: n a {r,t) = J 7^^/^(P> r > *)> 



(2tt/i) 3 - 

total density : n tot (r,t) = ^~]n a (r, t), 



velocity: n a (r, t)v a (r, t) 



dp p 



(2vr^) 3 m 

pressure tensor: P^(r, t) = ^^(r, t) 



(Al) 
(A2) 

(A3) 



E 



m 



dp 



— -v a ^(r,t) — - <u CT) „(r,t) f a (p,r,t), 



m 



.m 



energy density : E(r, t) = ^ E <r( r , 

a 

dp 1 



=E 



(2vr/i) 3 2m 
heat current : Q(r, t) = ^^Q (T (r,t) 

dp 1 



[p-mv ff (r,t)f/ ff (p,r,t), 



(A4) 



(A5) 



E 



[p - mv ff (r,t)]' 



rate-of-strain tensor : D a ^ lJ (r,t) 



(2ttH) 3 2m 
dv a ^{r,t) dv at „(r,t) 



P 



m 



v CT (r,t) f a (p,r,t), 



dx. 



dx, 



(A6) 
(A7) 



Indexes /z and ^ are Cartesian components. 

We assume that local velocities of two components are the same t> M (r,t) = f| iAt (r,t) = 
W| iM (r, t). This means that a rate-of-strain tensor is the same for two components, and hence 
we define D^ij.t) = D a ^ u (r,t). With the above quantities, generalized hydrodynamic 
equations are given by 



d 

— n a {r, t) + V r • K(r, t)v(r, *)] = 0, 



(A8) 

^: P ^ (r ' t] ~ [ ^ T(r ' t)ni(r ' t)] ' (A9) 



mn tot (r,t) —+v u (r,t)— v^{r,t) -- 
^E(r, t) + V r Q(r, t) + V r [E(r, t)v(r, t)] + ^ D^(r, t)P^(r, t) = 0. 



(A10) 



fJ,U 
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These hydro dynamic equations are obtained by multiplying Eq. ([T]) by 1, p and p 2 and 
integrating over p. The collision integral in Eq. ([1]) vanishes owing to the conservation law. 

In the collision-dominated regime, the first approximation to the distribution function 
is the local equilibrium distribution / CT (p, r, t). In local equilibrium, the hydrodynamic 
quantities are given by 

1 



Q(r,t) = 0, 



where A(r, t) is the local thermal de Broglie wavelength: 



(All) 
(A12) 

(A13) 
(A14) 



A(r,t) = 



27Tfi 2 



mk B T(r, t) 



1/2 



(A15) 



J~n( z (r{ r ,t)) is the Fermi function [47| given by 

1 



dx- 



x 



n-l 



(A16) 



T (n) J exp (x)z a 1 + 1 ' 

where T(n) is the Gamma function. With the above quantities, hydrodynamic equations in 
local equilibrium are given by 



d 

— n CT (r, t) + V r [n a {r, t)v(r, t)] = 0, 



(A17) 



mn tot (r,t) 



d_ 

di 



v(r,t) • V r 



^(r,t) + |v r E(r.mr.t) =vfr.f)- Y v P(r.t) 



v(r, t) = -V r • P(r, t) - V [^ T (r, t)n i (r, t)] , (A18) 

(A19) 



In order to treat departure from local equilibrium, we introduce the following form of the 
distribution function: 



f a (p,r,t) = f a (p,r,t) + f a (p,r,t) l-/ CT (p,r,t) *„(r,t) 



(A20) 



Since the number of particle, the total momentum, and the total energy is conserved, the 
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following three constraints are imposed; 

dp 



£7 

5^ / (2nh) 3 



(27r ^ 3 /a(p,r,t) l-/ ff (p,r,t) * (r (p,r,t) = 



r/p p 2 /.(p,r,t) l-Z^r,*) tf ff (p,r,t) = 0. 



(A21) 
(A22) 

(A23) 



The local equilibrium distribution ([3]) satisfies the detail balance of the scattering 
1 - / CT (1)] [l - /_ ff (2)] /_„(3)/ ff (4) = / CT (1)/- CT (2) [l - /_ CT (3)] [l - 7,(4)] . With a use of 
this relation, the collision integral in the right hand side of the Boltzmann equation reduces 
to 



Z coll [/ CT (l)] =L a [* a (l)] 

_2irg 2 f dp 



h 



dp 3 



x 



(2nhy J {2.Hf J dP45(Pl + P2 " P3 " P4)5 + £ - t - t 
1 - / ff (l)] [l - /_ CT (2)] /_ CT (3)/ CT (4) [^(4) + *_ ff (3) - #_ CT (2) - . 

(A24) 



It is useful to introduce dimensionless momentum variable £(r, £) = u(r, i)y m/[2/cBT(r, t)] 
where mu(r, i) = p — mv(r, £), and to introduce the dimensionless collision operator 

x 1 - £(1)] [l - /_ CT (2)] /_ CT (3)jV(4) [^(4) + *_ CT (3) - vP_ CT (2) - . 

(A25) 

The collision integral ^[^(l)] is then reduced to ^[^/^(l)] = L' a [^(l)] /C(r,t), where 
the coefficient C(r,t) is defined as C(r,t) = ti 3 n 3 / ^4a 2 m k^T^r^t) j. 

Here, let us substitute the distribution function in the local equilibrium to the left hand 
side of the Boltzmann equation: 



- + ^--V r - W ff (r,t)-V p 
at m 



/a(p,r,t) 



I ( d p „ \ , s mu 2 (r,t) ( d p „ , ... 

+ - • V r z g (r,t) + - J;/ M — + ii ■ V r T(r,t) 



mu(r, t) / <9 p 



2A; B T 2 (r,t) m 



*bT (r, t) V* + m ' V 'J " (r ' *> + i)VrC, ' (r ' " ' U(r ' 4) 
l-/a(p,r,t) /^(p.r,*). 



(A26) 
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This equation can be written in a simpler form as shown below. Note that density and 
pressure satisfy the following equations: 



dn a (r, t) 3 n ff (r, t) dT(r, t) _ 7(T (r, t)k B T(r, t) dz c {r, t) 



dt 



2T(r,t) dt 
_ 3 n ff (r,t) 
2T(r,t) 



2«r(r,t) 



5f 



(A27) 



v(r, f) • V r f (r, f) - 7.(r,t)fc B T(r,t) _ ( ^ _ ^ 



z a (r,t) 



E dPg{r,t) 



n a {r,t) [V r -v(r,t)] 
^ Pa{r, t) 



(A28) 



dt 



V r T(r, t) ■ v(r, t) — V r * ff (r, t) ■ v(r, i) 



2T(r,t) 
Jp.M) [V r • v(r, *)] }► , 



z CT (r,t) 



aP CT (r,t) 5P a (r,t)df(r,t) , n ff (r, t)fc B T(r, t) <9z CT (r, t) 



9* 



2 T(r, t) dt 



Zo-(r,t) 



9* 



m ^ / \ 5-P(r(r,t)„ ~, v ^o-fr, t)k B T(r, t) „ 



(A29) 
(A30) 
(A31) 



where we define j a (r,t) = J r i/2(z a (r,t))/[k B T(r,t)A 3 (r,t)}. Here, we shall consider the 
following equation: 

h a (r,t) dh a {r,t) dP a (r,t) 



^ la {r,t) dt Ot 
Using Eqs. (1X271) - (IA301) . and flA32l . one obtains 

= _ V (r, t ) ■ V r T(r, f) - If (r, f) [V r • v(r, f)] . 
From Eqs. (1A27j) . (1A31h and (1A33[) . one also obtains the following equation: 



(A32) 



(A33) 



dz a (r,t) 
dt 



-v(r,t) ■ V r z a (r,t) 



z a (r,t) 



n a (r,t)k B T(r,t) 



v(r,i) 



V r P CT (r,t)-|5^V r T(r,t) 
2T(r,t) 

(A34) 



Using the above equations, we reduce the left hand side of Boltzmann equation to 

U{p,r,t) 
mu 2 (r,t) 5%(z ff (r,<))' 



— + --V r -VC/ ff (r,t)-V p 
ot m 



T(r,t) 

777. \ - ^ . . 



k B T(r,t) 



2k B f(r,t) 2F 3/2 (z a (r,t)) 
u^r, t)u u {r, t) - 5^u 2 {r, t) 



-^d a (r,t) ■ u(r,i) 



X 1 - fa(p,r,t) f a (p,T,t), 



(A35) 
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where we define do-(r, t) as 
1 



d CT (r,t) = 



/c B T(r,t)n tot n to t 



VP ff (r,t) 



+ Vt/ ff (r,t) 



VP_ CT (r,t) 



Vf/_ ff (r,t) 



(A36) 



In the population balanced gas, one finds d a (r,t) = 0. As a result, the left hand side of 
the Boltzmann equation under the local equilibrium in the population balanced Fermi gas 
is reduced to 



+ V r - VU a (r,t)- V p 
at m 



m 



f(v,t) 



/<t(p, r, t) 



2 5 ^ 5 /2(^(r,t)) 
2F 3/2 (z <r (r,t))\ 



+2^Pv(r,t) 



&(r,t)&(r,t)-<S^-£ 2 (r,t) 



l-f a (p,r,t) f a (p,r,t). (A37) 



We introduce an ansatz for the departure from the equilibrium, which is 



vMO =C(r,t) 



2k B f(r,t) 



1/2 



+2^£V(r,t) 



VT(r,t)-g(r,t) 
f(r,t) 

1 



MO 



^(r,^(r,t)-^^ 2 (r,t) 



5.(0 ■ 



(A38) 



This comes from a consideration that the solution must be a linear function of VT(r, t) 

and D^ u (r, t), based on Eq. (IA37I) . We substitute this into the collision integral on the right 
hand side of the Boltzmann equation. Comparing Eq. (1A37j) with this result, one obtains 
the following relations: 



2 5 ^ 5 / 2 (^(r,t)) 
2^ 3/2 (^(r,t))J 



l-/a(p,r,t) /„(p,r,t), 



(A39) 



£ M (r,t)£,(r,t)--^ 2 (r,*) 



^(r,^(r,i)--^ 2 (r,i) 



l-/a(p,r,t) /^(p,^*). 



(A40) 



THe ansatz (1A38D automatically satisfies two constraints (1A21j) and ()A23j) . For the 
constraint (IA22p to be satisfied, the function A a (£) should satisfy 



£ / de(r, t)U{t, r, t) \l - r, t)} £ 2 (r, t)4,(0 = 0. 



(A41) 
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Transport coefficients such as the thermal conductivity and the viscosity are obtained 
using the ansatz (IA38I) . The thermal conductivity k is defined by 



Q(r,*) = -«(r,f)VT(r,*)- 



(A42) 



From Eqs. (jMj) . flA20|) . and (jMHJ) , it is given by 

1 1/2 



re(r,*) 



24a 2 



2k B T(r,t) 



m 



dt?{r,t) l-/ ff (p,r,t) / ff (p,r,0Ar(O- (A43) 



The shear viscosity 77 is defined by 

P(r,t) = c^f%(r,t)-27/(r,t) 



D Miy (r,i)--TrL>(r,t)^ i/ 



(A44) 



From Eqs. flA"20|) . and flA"38|) . it is given by 



7?(r,f) 



60a 2 



2k B f(r,t) 



rn 



1/2 



£ / ^£ 4 M) [l -; CT (p,r,t)l f a (p,r,t)B a (£). (A45) 



Note that the second viscosity (the bulk viscosity) is absent. In more general, the second 
viscosity vanishes in the normal gas interacting with the s-wave scattering, because a uniform 
compression at a steady rate changes the thermodynamic equilibrium into a new one (see 
the second of Ref. fl). 

Using these quantities, we reduce hydrodynamic equations for the velocity and the energy 
density to 



mn tot {r,t) 



d 



d_ 

di~ 
_d_ 

dx,. 



v(r,t) • V 



2 V {r,t) 



D^(r,t) - -TrD(r,t)5^ 



-E(r, t) + V [E(r, t)v(r, t)} + [V • v(r, t)} P(r, t) 



V[K(r,t)VT{r,t)] + 2 V (r,t)J2 



D^t) - -TrD{r,t)S^ 



(A46) 



(A47) 



Based on Refs. 



44 



46], we shall take the function A a (£) as 



Mt) = A 



5^5/2^) 



2^3/ 2 (^). 



(A48) 
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Note that this satisfies a constraint in Eq. (jA41B . Multiplying Eq. (1A39|) by £[£ 2 
|^ : 5/2(^(t)/^ : 3/2( ;Z (t)], and integrating over £, one obtains the coefficient A as follows: 



.4 



t2 _ 5 ( 2 <t) 
5 2^ 3/2 ( 2ct ) 



[1 



EM 



/•2 _ 5 ^5/2 (fcr) \ £ 
? 2^ 3 /2(« CT ) ' ? 



-2 5 5/2(^0-) 



15 7T 3 / 2 

T x7 



E 



7^7/2 Oa) 



where Xa is defined as 

2a = E / ^ 



.2 _ s fyiiz*, , 

s 2^ 3 /2(^) ^ 



5 Tzjiizo, 



5 ^1/2^°) 

(A49) 



(A50) 



On the other hand, based on Refs. 



44 



461 ]. we shall take the function B a (£) as -B<t(£) = 



£>. Integrating over £ and summing over a, z/ and /x, after multiplying Eq. flA40p by 
(£1^ ~ |<W£ 2 )' one obtains the coefficient S as follows: 



EEJ^ fc(r, t)£,(r, t) - §« 2 (r, t)] 2 /.(p, r, t) 1 - /.(p, r, f) 



5 



(T /if 



EEM fo(r,*)^(r,*) - |« 2 (r,t)] V 9 [fo(r, t)^(r, f) - f« 2 M))] 

CT /if 



(A51) 



where 

X B = ^ / d£ 



£ M (r,t)&(r,t)-^£ 2 (r,t) 



^(r,t)^(r,t)-i^ 2 (r,t) 



The collision integral satisfies the hermitian property 

E / ^(0^(0] = E / wo^iMa 



(A52) 



(A53) 



The collision integral in Eq. (IA25I) also satisfies -£4[£i] = 0' ano - = owing to the 

conservation of the momentum and the energy; then, Xa and Xb are reduced to simpler 
formulae given by 



j a = E / d ^ 2 K [en] , j b = E / d zuA \ua ■ 



(A54) 



We shall introduce new variables given by £ x = (£ + £') / a/2, £ 2 = (£ — £') / \/2, £ 3 = 
(£0 + £") / y^, and £ 4 = (£0 — £") / v^2- Note that these variables satisfy relations £ = £q 
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and = \$"\, because of the conservation of the momentum of the center of mass and of 
the energy in the relative motion. Here, we shall define a function F a (£o,£',y',y") given by 



F a (^,e,y',y") = Mi)f-a(2) W-,(3) [1-7,(4) 



(A55) 



[l + z a exp (-£?)] [1 + z_ a exp (-£ 2 2 )] [1 + z_ a exp (-££)] [1 + ^ exp (-£)] ' 

(A56) 

where £? = + ^ofV + £' 2 )A & = - 2£ £V + t 2 )A £f = (£o + 2^Y + £' 2 )/2, and 
£| = (£2 _ + £' 2 )/2. The function F a (£ , f, 2/', y") satisfies the following relations: 



E F <^o, -y, = E F ^ £'> ^. -y")- 



(A57) 
(A58) 
(A59) 



The integrals in ( 1A54I) in terms of the new variables are thus rewritten as 

/*oo /*oo r\ /*1 

lA=V2n* «o 4 / #W dy' dy"J2mo,e;y',y"W 2 + y" 2 -2y' 2 y" 2 ) 
Jo Jo J-i J-i 



(A60) 



-i poo POO pi pi 

is =^ 3 y o « 2 j o d^r J h J ^ d V " e e-, y', 10a + </ 2 + y" 2 - s y 



(A61) 



For simplicity, we again introduce variables £o = ^/2rj cos and £' = y/2rj sin 0. The 
condition of the population balanced gas is given by fa — p,±. In this condition, the function 
F a (rj, 0, y', y") can be reduced to 

,11 1 



4 cosh (rj — (3p, to t/2) + cosh (7/2/ sin 20) cosh (77 — /3/j tot /2) + cosh (777/" sin 20) 

(A62) 
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where /i to t = A*t + AH- As a result, integrals X4 and Z# are reduced to 



11/2 



J4 =2V d / drm 

'0 -/ll 



' cos 4 sin 7 1 



J^dy' J^dy"(y' 2 + y" 2 -2y' 2 y" 2 ) 



1 



x 



cosh (77 — /5/i to t/2) + cosh (r/?/ sin 20) cosh (77 — (3j2 tot /2) + cosh (777/" sin 20) ' 

(A63) 

Jb=2 3 7T 3 / dr/r/ 9 / 2 / d0 cos 2 sin 7 / <V / + y' 2 + t/" 2 - 37/ ,2 7/" 2 ) 

Jo io J-l J-l 

1 1 



X 



cosh (77 — Pntot/2) + cosh (777/' sin 20) cosh (77 — (3jl tot /2) + cosh (777/" sin 20) 

(A64) 

Before closing this section, we summarize that the thermal conductivity k and the vis- 
cosity 77 are given by 



k = _ kv f 2fc B T(r,t) ^ V2 75^3 J_ I J2 



m 



128 J A 



7 5 ^ 5 2 /2 (^ 
2^ 2 ^ " 2^M^ 



1/2 



m / 2£; B T(r,t) \ 5 3 1 
V ~ ~ [ m J 32^ 2b 



(A65) 
(A66) 



where Eqs. flA43|) . flA45|) . (IA491) and flAoll are used. 



APPENDIX B: RELAXATION TIME 



A purpose of the present appendix is to derive relaxation times using results obtained in 
the previous section. Let us consider the solution in the collisional hydrodynamic regime. 
In this regime, the departure from local equilibrium on the left hand side of the linearized 
Boltzmann equation (Tl9l) are neglected. Solving it for 5u cr (p, q, o>), one obtains 



(Bl) 



<5z/ CT (p, q, lo) = it^(uj- ^—^j [a CT (q, u) + b(q, to) ■ p + c(q, u)p 2 ] 

£ [a_ CT (q, lu)W^ + c(q, uj)W^ 2 ] } . 



A closed set of equations for a CT (q, to), b(q, to) and c(q, u) can be obtained from Eq. (JTH 
by multiplying Eq. (fl9l) by 1, p, p 2 and integrating over p. The zeroth moment yields 



uja a [q,LO) 



W*2 
W- 



b(q,^) • q 
3m 



(B2) 
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The second moment yields 



wc(q, to) — 



b( q^) • q 

3m 



(W T)4 + W i>4 ) = 0. (B3) 



From Eqs. (jB2J) and (JB3J), relations given by 

b(q,w)-q <Jv(q,w)-q 

uc(q, (jj) = = , 

3m 3m 

and ao- (q, oS) = are obtained; therefore the departure from local equilibrium in Eq. flBT 



(B4) 



is reduced into 



x ( \ ■ p ' q « 



V 



w n 



W r 



(7,2 



59(q, lo) + ir 



<Jv(q,w) q 2 p q , 

-p + — dv(q,w) • P 



3m 



(B5) 



From Eq. (|A42j) . the heat current in the Fourier representation is given by Q(q, to) 
-/dq<5T(q, u). From Eqs. (]A6j) and (1B5I) . the thermal conductivity k is thus obtained as 

for 



K = — 



12m 4 



" (7,4 



(B6) 



The Fourier representation of the rate-of-strain tensor in Eq. (1A7j) is given by D M „(q, u; 
i [^^(q, u;) + <?yW M (q, a;)] /2, and hence the pressure tensor is P fU/ (q, a;) = <5 Mj! ,P(q, u; 
277 [P Mi ,(q, u;) — 5 |t)I/ Tr£)(q, a»)/3]. From Eq. (1B5I) . the viscosity 77 is obtained as 

rW aA 



V 



E 



15m 



2 ' 



(B7) 



Comparing Eq. (1A65|) with Eq. (1B6|) . we obtain the relaxation time associated with the 
thermal conductivity (denoted by r K ) given by 



15 



7T 



13/2 H 



g 2 m 3 (kBT) 2 2a 



--It 



7 5 -P 5 2 / 2 (^) 
2^ ( ^ " 2 7^) 



4vr 5 ^ 7 



g 2 m 3 (kBT) 



-A. 



(B8) 



Comparing Eqs. (1A66P and (lB7p . on the other hand, we obtain the relaxation time associated 
with viscosity (denoted by r v ) given by 

47T 5 H 7 



g 2 m 3 (kBT) 2 Xb 
We have used the following relation: 



g 2 m 3 (kBT) 



:B. 



(B9) 
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Equations flA65|) and (IA661) are written in terms of the temperature and the fugacity in 
local equilibrium. We note that these quantities should be taken as equilibrium values in 
the expressions for the relaxation times. 

On the other hand, the mean-collision time r co n is defined by 



N tot 1 



(BH) 



V r coll J (2vr7i)3 

lug 1 f dpi f dp 2 f dp 



dp^Pi + p 2 - P3 - P4) 



H J {2-nkf J (2ixhf J {2nhf 

Following the procedures analogous to those deriving Eqs. (IA63I) and (IA64[) . we reduce r co n 
to 

N tot 1 2 m 9 / 2 (fc B T) 7 / 2 /"» 5/2 r'\ . 3 f 1 f 1 

1 1 

X 



cosh (77 — (3fj, tot /2) + cosh (r/?/' sin 20) cosh (77 — (3fi tot /2) + cosh (777/" sin 20) 

(B13) 

In Fig. [51 viscous and thermal relaxation rates are plotted. The mean-collision rate is also 
shown. The coupling constant a = gN tot /VeF = 1 is used, where Ef is the Fermi energy. 
Behavior of the viscous relaxation time is different from that of the mean-collision time, as 
noted in Ref. [ijj]. The viscous relaxation rate is severalfold bigger than the mean-collision 
rate, and is effective in the hydrodynamic regime compared with other relaxation rates. In 
the low temperature regime, although the thermal relaxation rate is bigger than the viscous 
one, the difference is very small. 

In summary, the viscous relaxation rate is the most important in the high temperature 
regime. We apply this viscous relaxation rate to the relaxation time in the moment method, 
because the density oscillation is the most strongly coupled with the viscous relaxation, and 
this relaxation rate is dominant in the high temperature regime. 



APPENDIX C: RANDOM PHASE APPROXIMATION 

We solve the linearized Boltzmann equation in the collisionless limit using the random 
phase approximation. For this purpose, we add a small perturbation U a (q,u) to evaluate 
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0.7 




k B T/£ F 



FIG. 5: Relaxation rates 1/r versus temperature. Viscous and thermal conductivity relaxation 
rates are shown with solid and dashed lines, respectively. The mean-collision rate is plotted with 
dotted line. A coupling constant gN tot /V is assumed to be the Fermi energy: a = gN tot /VeF = 1. 



{(^ ~ ^( q ' P '^ + ^^'^ + 9^ a (q,uj)]^ = 0, (CI) 



the density response function. The linearized Boltzmann equation becomes 
del 

where we neglect the collision integral on the right hand side. The fluctuation around static 
equilibrium is thus given by 

v a (r, p, uj) = ^p^ P q [^(q> P' w ) + 9 s n- a {(l, uj)} ■ (C2) 

m 

Since the density fluctuation can be written as 

5n ff (q,u) = J ^^(q.ftw), (C3) 

the density fluctuation in terms of a response function x^i^ 00 ) * s gi yen by 5n a (c{,u)) = 
X°(q, uj) [U a (q, uj) + g6n- a (q, uj)], where the density response function x° (q, to) is defined as 

o/ \ f d P d fa 1 P-q f n ^ 

We assume that the density perturbation is the same for the two components: U(q, uj) = 
U a (q,uj) = C/_ CT (q, uj). The density fluctuation is then reduced to Sn a (q, uj) = 
x(q, uj)U (q, uj) , where the response function is given by 

i - gxl (q> "JX-rfa OJ) 
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Zero of the denominator of a response function gives the frequency of the collective mode. 
In this case, the dispersion relation of the zero sound is obtained from the following equation: 

l-g X \c l ,uj) = 0, (C6) 

where we assumed the population balanced gas and used the relation u ) = w ) = 

Note that, this equation ( 1C6I) at T = reproduce the dispersion relation ( 1561) . Solution 
of the linearized Boltzmann equation involves the denominator uj — p • q/m, as seen in Eq 
( 1C4|) . This means that excitations within the linearized Boltzmann equation can reproduce 
only the phonon regime: Q oc q. 
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